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1 Introduction 



In this paper we consider the mass-shell and the Galilean limits of equilibrium rel- 
ativistic distribution for indistinguishable events studied in a previous work Q. In 
that work we studied an identical many-body system within the framework of a man- 
ifestly covariant relativistic statistical mechanics discussed in a series of papers [2]- [4]. 
In this framework, for an iV-body system, the N events generating the N particle 
world lines are considered as the fundamental dynamical objects of the theory; they 
are characterized by positions g M = (c£, q) and energy-momenta p M = (E/c,p) in an 
8iV-dimensional phase space. Their motion is parametrized by a continuous Poincare- 
invariant parameter called the historical time. Such a system of N events is described 
by generalized Boltzmann equation [5], whose equilibrium solution gives the distri- 
bution functions (for both bosonic and fermionic events treated simultaneously in 
one expression) coinciding with the corresponding grand canonical distributions ob- 
tained in ref. [2] for the static Gibbs ensembles. Upon integration of this distribution 
function over angular and hyperangular variables, one obtains the relativistic mass 
distribution [1]. We found expressions for the pressure and the energy density in such 
a system and obtained the relativistic equation of state. Now we turn to a mass-shell 
form of that equilibrium distribution and to its Galilean (nonrelativistic) limit. 

The Galilean limit of a manifestly covariant relativistic statistical mechanics was 
considered in refs. 0,0, by taking c — > oo (compared to all other velocities). In this 
limit the relativistic relation between the energy E and the mass m 

E 2 = m 2 + p 2 

transforms to 

P 2 

E = m + — . (1) 
2m 

If we require in addition [6], [7], that the quantity 

n = c 2 (m - M) (2) 

may take any value, however, finite, as c —>■ oo, then m = M(l + 0(l/c 2 )) 

(i.e., the "mass-shell limit"), and the relation between E and m takes on the form [7] 

E = m + Wr < 3) 

where the Galilean mass M coincides with the particle's intrinsic parameter. 

In case of an equilibrium relativistic ensemble of distinguishable events such a 
transformation of the relativistic relation between E and m gives rise to the usual 
nonrelativistic Maxwell-Boltzmann distribution of p 2 /2M [7]. 

In this work, we show that in the mass-shell limit, the equilibrium relativistic 
distribution for indistinguishable events used in ref. [1] approaches the distribution 



2 



found by Jiittner and Synge M within the framework of an on-mass-shell relativistic 
kinetic theory. We also study the c — > oo limit of the on-mass-shell theory, and show 
that it goes smoothly to the Galilean form. 



2 Equilibrium relativistic distribution close to mass 
shell 

Consider the equilibrium relativistic distribution for indistinguishable events used in 
ref. [1] (we use the metric g^ u = (— , +, +, +) and q = g M , p = and assume no 
degeneracy; in the case of degeneracy all the corresponding formulas throughout the 
paper should be multiplied by degeneracy factor), 

e A(q)(p+p c f+B(q) 

f (q,p) = C( q ) lTeMq)ip+pc)2+B{qy A(q)>0, (4) 
which is normalized as 

d 4 pfo(q,p) = n(q), (5) 



where n(q) is the total number of events per unit space-time volume in the neighbor- 
hood of the point q. 

It follows from the relations [7] (we suppress c for the present consideration) 

i] = m — M, -A <r]< A 

that 

M-A<m<M + A. (6) 

Since A may take any value as small as one wishes,Q but not zero, i.e., the variation 
in mass of the particles of the ensemble may be very small, we can take the value of 
p 2 = —m 2 restricted to a small neighborhood of a fixed value — M 2 . This permits us 
to write (4) as || 

e -A(M 2 +ml)+B ^.Ap^vcy. 
/OWJ-P) — C(<Z) ^ e -A{MZ+ml)+B e 2Ap»p cl , ' ^) 

Introducing hyperbolic variables [3] and performing integration we obtain from 
(5) and (7) the normalization relation 

Air AM 2 00 f-M N i n + 1 

n ( q ) = C(q) ~~ V ^ e~ nA ^ M 2 +^)+^^ (2nAMm c ) , (8) 

Am c ^ n 



1 It corresponds to the approximate mass shell condition | m — M |< A. 
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where K\ is a Bessel function of the third kind, defined by 



K v (z) = ^e^H^(iz). (9) 

The average values of p^, p^p v ', etc., can be obtained from the corresponding 
relations of ref . [1] , using the formula ITT 



TL 

K' n (z) = -K n+1 (z) + -K n (z). (10) 

Z 

Thus 

^ _ „ M E^ = i{(±l) n+1 e- nA(M2+m ' )+nB K 2 (2nAMm c )/n} 
{P }q ~ Pc m c E"=i{(±l) n+1 e- nA ( M2+m ^ +nB ^i(2ny4Mm c )/n}' { } 

M ^ M En=i{(^) n+1 e- nA(M2+m ' )+nB K 2 (2nAMm c )/n 2 } 

{P V )q - 9 2Amc J2^ =i {^n+l e -nA(M^rnl)+nB K ^ 2rl AMm c )/n} 

, M 2 gM EgL 1 {(±l) n+1 e-^ (M2+m - )+ " i? ^3(2nAMm c )/n} 

E~ = i{(±l) n+1 e-" i(M2+m ' )+,l ^i(2nMm c )/n}' 1 j 



Identifying, as in [1][|, 

2Am 



1 



M 



m. 



B = ]£t (" + t)- (13) 

where T, /i and are absolute temperature, chemical and mass potentials, respec- 
tively, we have 



n/j, 



E^-i{(±l) n+1 e^^ 2 (^)/n 2 } 

(pV) = g^MkgT — , kBT 

Y^ =l {{^Y^BT Kl {^)/n} 

I 

TbU 

^^ ^{(iire^^i^l/n} ^ (i5) 
K Y% = A^Y +1 ^K x {^)/n} 



2 The first relation implies that in thermal equilibrium Am c is independent of q. 
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where \x = /i — (N being the total number of events) is the "reduced" chemical 
potential approaching /iasiV-^ooorT-^0 (since in the latter case fix ^ [1]). 

As in previous works [1],[3],[4], to obtain the local energy density we make a 
Lorentz transformation to the rest frame of the local average motion. According to 
(14), the relative velocity of the new frame is 

Pc 
U = . 



The rest frame energy is then 



{E)= (E)-u- P 



2 



so that 



E~i{(±l) n+1 e^^ 2 (f^)/n} 
(E) = M , VkBT " . (16) 



E^ =1 {(±l)" +1 e^^(0)/n} 

To obtain the pressure and the energy density in our ensemble, as in previous 
works [1],[3],[4], we study the particle energy- momentum tensor defined by the R 4 
density, 



T^(q) = T,Jdr P -^5\ q - qi (r)). (17) 

i 

Using the result of [3] 

(T^(q)) q = T AV Jd 4 pfo(q,p)^f (18) 
and the expression (15) for (p^p u ) q , we obtain 



{T ^ q))q = T Av n( g) 



M 



' M Y,Zi{(±l) n+1 e- nAiM2+m * )+nB K 2 (2nAMm c )/n 2 } 
9 2Am c Y.n=i{( ±l ) n+le ~ nA{M2+ml)+nBK i( 2nAMm c)/n} 



PKE™ = i{(±l) n+1 e" ni(M2+m ' )+nB ^3(2nAMm c )/fj} 



m 2 



E~ = i{(±l) ri+1 e- nA ( M2+m ') +nB Ki(2nAMm c )/n} 



(19) 



In this expression T&y is the average passage interval in r for the events which pass 
through a small (typical) four-volume AV in the neighborhood of the point q of R A . 
The formula for the stress-energy tensor of a perfect fluid has the form [3] 



where p is the pressure and p is the density of energy at q. 
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According to (14), 



!~{p x )M q mc 

hence 



(21) 



_ T Av n(q) J2^ =1 {(±l) n+1 e- nA ^ M2+m ^ +nB K 2 (2nAMm c )/n 2 } 
P ~ 2Am c E"=i{(±l) n+1 e" ni(M2+m?)+n ^i(2nMm c )/n} ^ ^ 



and 



E^ 1 {(±l)" +1 e-- 4 ( M2+ ^) + " B K 3 (2nAMm c )/n} 
p + p _ ^Ayn( g jM^ {(±i)n+le _ nA(M2+mi)+nBXi(2 ^ Mroc)/ ^. ^dj 

To interpret these results, as in [1],[3],[4], we calculate the average (conserved) particle 
four-current having the microscopic form 

nq)=T,Jjj^-^))dr. (24) 

i 

Using the result of [3] 

(nq)) q = T AV Jd 4 p^f (q, P ) (25) 

and expression (14) for (p^) q , we obtain 

HHaW -T n (g\ P " Sr =1 {(±l)" +1 e- i(M2+m - )+nJj ^(2nM mc )/ n } 

{ {q))q AV Wj m c E^ = i{(±l) n+1 e-^( M2 +^)+" B ^i(2nAMm c )/n}- 1 } 

In the local rest frame p^. = (m c , 0), 

E- =1 {(±l)" +1 e-^ M2+ ^) + "^ 2 (2nAMm c )/n} 
W - J W^ E ~ 1 {(±l)«+i e -^(^ 2 +^)H-^ J ft: 1 (2nAMm c )/n}- 1 } 

Defining the density of particles per unit space volume as 

N (q) = (J°(q)) q , (28) 

we obtain the relativistic equation of state 

N En=i{( ±1 ) n+1 e~ nA(M2+m2c)+nB K 2 {2nAMm c )/n 2 } 



p = 



2Am c J2^ =1 {(±l) n + 1 e~ nA ( M2 + m2 c)+ nB K 2 (2nAMm c )/n} 
Xn=A(±l) n+1 e^K 2 (^)/n 2 } 



N k B T — — — -. (29) 

' nM \ 
-k B T> 



E~=i{(±l) n+1 eW 2 (^)/n} 
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It follows from (22), (23) that 



n/j, 



p = M E^ =1 {(±l)" +1 e^^3(|^)M _ 1 
P En =1 {(±^ n+1 e^K 2 (^)/n^} 
Let us introduce T through the relationP] 

p=(T-l)p. (31) 

Then 

r = i + —, . (32) 

Er= 1 « ±1 ) n+le ^ r ^( 1 ^)/n} 



M - -,.■■-.1.' ~*K B J. " J -J^ 



Using the asymptotic formulas [12 



^-e z , z — > oo, 
Kv(*)~{ L „ (33) 



|r(i/)(|)", ^o, 



one obtains 



1 T — > 

r= 4 i ' (34) 
\ 3, T^oo. 

Since in thermal equilibrium there is no dependence on q, instead of fo(q,p) one can 
use equilibrium relativistic distribution function [5] 

/o(P) = 7Z ~ ~ , (35) 

normalized as^ 



(25T)' 

where n = is the event number density. In this case one obtains, instead of (8), 



~! n \k B T J 



n 



and then, from (22), (23), (27), 



3 This relation is helpful when studying an adiabatic equation of state |n 
4 Eq. (36) is formally equivalent to (4), (5), with C(q) — . 
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n=l 

Using now the relation [TH 



2n 

K n+1 (z) - tf„_i(z) = — tf n (z), (41) 

z 



one has from (38), (39), 



(T AV A)M 3 k B T ~ (±l) n+1 



The three relations (38), (40), (42) coincide with the corresponding relations of an on- 
mass-shell relativistic kinetic theory (see Appendix A), except for the additional factor 
(TAyA)/7T. Therefore, (both the on-shell and the off-shell in the limit m 2 = M 2 ) 
approaches give the same results, if one takes T A y/\ = ir. Since 2A = 5m is the 
width of the mass deviation from its sharp value m 2 = M 2 , the latter relation can be 
rewritten as 

T AV 5m = 2tt. (43) 

The presence of the factor {T AV A)/rr in the expressions (38), (40) and (42) shows 
that, as 5m — > 0, the corresponding quantities would go to zero, unless T AV —> oo. 
Preservation of finite values of p, p, N in this singular limit implies that the presence 
of events in the ensemble with very sharp mass requires the extension of the set of 
events over a very wide range of r, and therefore, t, since they are related through 



the equation of motion [15 



dt E 



hence 



IE) E~-i{(±l) n+1 e^ f ^ 2 (^)M 
St = (dt) = (dr) [ -f = T Av ^ n - lU Vk ° T)/ \ (44) 

where we used the relation (dr) = T A y and the formula (16). In the Galilean limit 
c — > oo the argument of the i^-functions goes to infinity^, so that one obtains, using 
the asymptotic formula (33) for z — > oo, 

(dt) = (dr), (45) 

in agreement with [6], [7]. 

5 This argument is actually equal to n ^ T ■ 
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If 5m is not so small, the range of r, and therefore, t, need not be very large. In 
this case, for each value of the historical time parameter r, the set of events making 
up the ensemble may be concentrated in a range of t which is not too large. The entire 
ensemble then moves in phase space with the development of r, and its stationary 
form corresponds to equilibrium. For 5m — > the ensemble fills a longer tube in 
phase space, and its displacement with r is not as important. It is in this limit alone 
that a direct comparison with the usual (on-mass-shell) relativistic ensemble can be 
made, since in this case the system possesses a stationarity in space-time but not a 
non-trivial evolution in r, as was remarked in ref. [4]. 



3 Limiting cases 



Now we wish to consider the expressions for N , p and p in the limiting cases M — > 0, 
T — ► oo, and T — ► 0, and to compare them with the corresponding formulas of an 
on-mass-shell relativistic kinetic theory, 
(i) M -> 0. 

We remark that the dynamical evolution of the system contains the quantity 
^jj- [15]. In the limit M — > oo, this quantity is undefined; however, if we first 
go to the mass-shell limit, for which there is no non-trivial evolution in r, then this 
quantity becomes numerical- valued (— M/2) and the limit M — > is well-defined. We 
therefore consider, for this limit, the on-shell results (38), (40), (42), obtained from the 
distribution function (35). 

Using the asymptotic formula (33) for z — > 0, one obtains the relations 



No = ± (T AV A)(k B T? Lk(±e 4 T . 



7T" 



P 



(T AV A)(k B TY 



7T" 



p = ±3 V , ; Li 4 (±e k BT] 



7T 



(46) 

(47) 
(48) 



(here Li v \z) = Y^k=\ fir is the so-called polylogarithm |]I6]|), which coincide with the 
corresponding relations of an on-shell theory, if T A yA = it (see Appendix B). 
(ii) T -> oo. 

This limit is mathematically equivalent to the limit M — > 0; physically it means 
that the particles become distinguishable and possess the relativistic distribution 
found by Synge |T7] . Thus, one obtains the same expressions without summation on 
k : 

(T AV A)(k B T) 3 J_ . , 

Na = V „ ; e fc a T , (49) 



7T" 
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V 



p 



(T Ay A)(^T) 4 



TV 



(T AV A)(k B TY 



g k B 1 



71 



(50) 
(51) 



g. 



in agreement with the textbook results, if T A yA 
(hi) T -> 0. 

We study here the case in which the theory is on-mass-shell. If, in the limit 
T — > we do not insist on the mass-shell constraint, then the mass distribution 
remains essentially relativistic. This case was treated for distinguishable events in 
ref. [7]. The indistinguishable case is treated in ref. [13]. 

Using the asymptotic formula (33) for z — ► oo, one obtains 



N 



T AV A (Mk B T* 



3/2 



7Y 



V 



T AV A (Mk B T 



2tt 



3/2 



Li 3/2 (±e k i 



-M 

7T' 



T AV A (Mk B r 



71 



3/2 



27T 



« P, 



71 



2tt 



k B TLi 5/2 (±e "b 

-M 

MLi 3/2 (±e "bt ) + 3k B TLi 5/2 (±e "b 



(52) 



(53) 



u —M 
' W s 



-- MN + 3p, 
(54) 

the formulas of a nonrelativistic statistical mechanics of indistinguishable particlesQQ, 
if one takes T A yA = ti and /i — M = /x — M = /i nr , the nonrelativistic chemical 
potential [H5| . 

As was remarked in [7], the low-temperature limit does not necessarily coincide 
with the nonrelativistic limit of a theory (for example, very long wave length radiation 
in Maxwell's relativistic theory does not coincide with the static Coulomb limit; it is 
still radiation). We see, that the procedure of taking the limit T — > on mass shell 
leads to a deformation of the symmetry group from Poincare to Galilean invariance, 
which means that this procedure actually coincides with the Galilean limit c — » oo 
PH (formally equivalent to the limit T — > 0, since in both cases the argument of the 
^-functions in the formulas (38), (40) and (42), z = , goes to infinity). We shall 
return to a more detailed study of this point in a future research. 



3 Indeed, using the formula (B.l) of Appendix B, one has (e 



2m > 



N 



d 3 p 



1 



,3/2 



± 



mksT 



3/2 



Li 3/2 (±e k Bi-), 



the normalization relation of a nonrelativistic statistical mechanics of indistinguishable particles, 



etc. 



20 



7 See the corresponding formulas of a nonrelativistic statistical mechanics of distinguishable par- 
ticles in ref. [18], p. 62, formulas (3.55). 
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4 Nonrelativistic distribution 

2 

In conclusion we shall show that in the Galilean limit the variable E — m = ^ has 
the usual nonrelativistic distribution for indistinguishable particles. 

We start with the normalization relation for the sharp-mass form m 2 = M 2 of 
initial equilibrium relativistic distribution (4), 

,. e -AM 2 -Aml+2App c +B 

n(q) = C(q) J d 4 p - - e _ AM ^ Ami+2ApPc+B ■ (55) 
This integral written in the local rest frame takes the form 

/oo r 00 e ' 

dEd 3 pY / (±l) n+1 e- 2nAmcE e-^ AM2+Am "- B) = / dEd 3 pJ2(±l) n+1 e-^e^, 

n=l •* n=l 

(56) 

or, in view of (2), (3) and the relation — A < rj < A, 

V / dme k s T d 3 pe fc s TlaM Mj . (57) 

n=l JM-A J 

Integration on m gives 

fM+A _rvm_ nM 

/ dm e "b t = 2Ae fe s T ; (58) 
Jm-a 



therefore 



1 



/ p ; f 
d 3 p E (±ir +1 e-^ ( ^ +M) = 2 AC(q) / d 3 p — r — , 
n=l J e (jM-v)/ k BT -j- 1 

(59) 

which is the usual (normalized) nonrelativistic Bose-Einstein/Fermi-Dirac distribu- 

2 

tion (we call e = ^) 

pl/2 

/( e ) = D , (60) 

ex P (f^)Tl' 



3 W ^3/2^ (±ir" 



^ = rg(^E^". (6i) 



In this formulas 



fi = li - M = f,-M[l + ^j (62) 

is the "reduced" nonrelativistic chemical potential approaching, as N — > oo or T — > 0, 
the value (/x — M) [1],[19]. 
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5 Concluding remarks 



We have considered the limiting mass-shell form of equilibrium relativistic distribution 
for indistinguishable events, studied previously in [1], and its Galilean limit. We have 
shown that the mass-shell forms of the basic thermodynamic quantities coincide with 
the corresponding expressions of a usual on-shell relativistic statistical mechanics, 
providing the important relation (43) for a width of the mass deviation over its sharp 
value. This relation implies that such a mass deviation in the system of events can 
be infinitely small. In this case the system of relativistic (mass shell) events becomes 
a stationary distribution in space-time, representing in this way the ensemble of the 
relativistic on-shell particles. 

In ensembles, in which the appreciable mass fluctuations can occur, the relativis- 
tic distributions (the equilibrium solutions of the generalized Boltzmann equation) 
represent the corresponding equilibrium relativistic distributions of mass, considered 
previously in refs. [1],[4], for identical and nonidentical systems, respectively. 

We have seen that the characteristic limiting cases give results in agreement with 
the corresponding limits of the usual theory; the low-temperature limit on mass shell 
in fact corresponds to the Galilean limit c — > oo, in which the statistical mechanics 
of indistinguishable events goes over to the usual nonrelativistic statistical mechanics 
of indistinguishable particles. 

For ensembles with mass fluctuations the framework of a manifestly covariant 
relativistic statistical mechanics provides corrections to results obtained within the 
usual on-shell theory. Physical consequences of the mass fluctuations for relativistic 
systems are considered in refs. [13l,[p^|,p5|, taking account anti-events, i.e., the 



events having the opposite sign of the mass potential \ik- In ref. [13] an adiabatic 
equation of state, p oc N^ 5 , is obtained for the system of degenerate off-shell fermions 
and possible implications in astrophysics are discussed. In [24] statistical mechanics of 
the bosonic event-anti-event system is considered. For such a system, at some critical 
temperature a special type of Bose-Einstein condensation sets in, which provides the 
events making up the ensemble a definite mass and represents, in this way, a phase 



transition to the usual on-shell sector. In ref. |25j possible applications of the off-shell 
theory in hadronic physics are discussed. The equation of state p oc T 6 , obtained in 
p5| , corresponds to the "realistic" equation of state, proposed by Shuryak |26| for hot 
hadronic matter, which is in good agreement with experiment for the temperature 
range 0.2-1.0 GeV. 
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Appendix A 



Consider the following relation, which represents the normalization condition for the 
relativistic distribution function 

/o(|p|,m) 



e (VP 2 +m 2 -/i)/fc B T -p 1 



(e = VP 2 + m2 ) : 



f d 3 p 1 1 /•«> e -( £ -^/ fcsT e v / e 2 - m 2 

~ J (2vr) 3 e (V?W-M)/feBT^7 ~^'m 6 1 =p e -( £ -M)/*BT • ( "J 

Expanding the denominator into power series and changing variables t = e/m, one 
obtains 

N = — - V(±l)" +1 e^ / dt t\fW^le~^ 1 . (A.2) 

27r n=l Jl 

Since [|3|] 

f* v^.-* = ^ (££) , (a.3) 

Ji nm \kb1 J 

one uses parametric differentiation with respect to (nm/ksT) and Eq. (10) of the 
main text to obtain 

r°° / nm. 8 rOO . nm, k rT [ 11171 \ 

/ dttVW^e-^ = --?— dtVW^le-^^^KJ^). (A.4) 
h d{^%) Ji nm \k R T J 



k B T> 



Thus, 



m 2 k B T ~ (±l) n+1 i, /nm\ , . . 

^ = -^J- E e*B" Jf a — . (A.5) 



JL 

2ir 2 ^ n " £ \k B T 



In the same way one can obtain [22 



/■ d 3 p p 2 1 _ m 2 (k B T) 2 ™ (±l) n+1 / nm 

P ~J (2vr) 3 3VP 2 + m 2 ^+^-/<)/*bT qz i ~ 2^ "i" ^ ^ 2 \A^T 

_ (A.6) 

f d 3 p VP 2 + ™ 2 m 3 fc B T ~ (±ir +1 / nm \ 

Zp+^r^l^— e" B T Kl — . (A.7) 



(2tt) 3 e (y/&+rt-Mk B T -p x 2vr 2 ^ n \fc B T 
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Appendix B 



First consider 

/•oo n u rln r°° p —a+b v 

r _^da_ = r daa „_e ^ = T,(±l) k+1 e kb f da 

Jo e a - b =F 1 Jo 1 =F e- a+b £r[ Jo 



oo 

a v e~ ka 



k=i K 



where Li n (z) = X^li fn: is the polylogarithm [16]. Then 

r d 3 p 1 4tt /• p> _ (k B T)* » 

7 (2vr) 3 edp|-M)/^ T ^ 1 (2tt) 3 7 e MW Tl * ^2 * j, 

(B.2) 

/ ^ 3 P P 2 1 4tt r p 3 tfo _ + (^) 4 rr+ ^ 

" ' (27r) 3 3|p|e(lpl-^)/ fc sT Tl 3(27r) 3 7 e^ p ~^/ kBT =F 1 tt 2 e s j, 

(B.3) 

r rf3 p IpI = _±_ f P'dp = ±3 (Ml Lu( ±e^) 
J (27r) 3 e(lpl^)/ fc sr T i ( 27r )3 j e (p-n)/k B T ^ ^ ^° ^2 J- 

(B.4) 
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